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Evaluating double and triple (?) boxes 
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Moscow 119992, Russia 

A brief review of recent results on analytical evaluation of double-box Feynman integrals is presented. First 
steps towards evaluation of massless on-shell triple-box Feynman integrals within dimensional regularization are 
described. The leading power asymptotic behaviour of the dimensionally regularized massless on-shell master 
planar triple-box diagram in the Regge limit t/s — > is evaluated. The evaluation of the unexpanded master 
planar triple box is outlined and explicit results for coefficients at , j=2,. . . ,6, are presented. 



1. Double boxes 

Fcynman diagrams with four external lines con- 
tribute to many important physical quantities. 
They are rather complicated mathematical ob- 
jects because they depend on many variables: in- 
ternal masses, Mandelstam variables and squares 
of external momenta. The most complicated two- 
loop diagrams are the planar double box and non- 
planar (crossed) double box. Almost all available 
analytical results correspond to the massless dia- 
grams. Ironically, the first result for the massless 
double boxes was obtained in the most compli- 
cated case, where all the fours external legs are 
off-shed, i.e. 7^ for aU i = 1,2,3,4. This is 
an elegant analytical result for the scalar planar 
master double box (and, moreover, for a general 
planar ladder diagram of this type), i.e. for all 
powers of propagators equal to one, obtained in 
[[^]. However, no other results for the pure off- 
shell double boxes (e.g. with dots on some lines 
and/or with some lines, other than rungs, con- 
tracted) have been derived up to now so that this 
result stays unique in the pure off-shell category. 

For massless double-box diagrams with at least 
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one leg on the mass shell, i.e. p| = 0, infrared 
and coUinear divergences appear, so that one in- 
troduces a regularization which is usually chosen 
to be dimensional [ ^ , with the space-time dimen- 
sion d as a regularization parameter. One hardly 
believes that a regularized double-box diagram 
can be analytically evaluated for the general value 
of the regularization parameter e = (4 — (i)/2, and 
the evaluation is usually performed in a Laurent 
expansion in e, typically, up to a finite part. 

The problem of the evaluation of Feynman in- 
tegrals associated with a given graph according to 
some Feynman rules is usually decomposed into 
two parts: reduction of general Feynman integrals 
of this class to so-called master integrals (which 
cannot be simplified further) and the evaluation 
of these master integrals. A standard tool to solve 
the first part of this problem is the method of 
integration by parts (IBP) [ |^ when one writes 
down identities obtained by putting to zero var- 
ious integrals of derivatives of the general inte- 
grand connected with the given graph and tries 
to solve a resulting system of equations to obtain 
recurrence relations that express Feynman inte- 
grals with general integer powers of the propaga- 
tors through the master integrals. 

The most complicated basic master planar and 
non-planar on-shell massless double-box diagram 
were calculated in [ ^, || by a method based 
on Fcynman parameters and Mellin-Barnes (MB) 
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representation (see [ g for details of the method) . 
It turns out that it is natural to consider non- 
planar double boxes as functions of the three 
Mandelstam variables s, t and u not necessarily 
restricted by the physical condition s + t + u = 
which does not simplify the result. 

Reduction procedures for the evaluation of gen- 
eral double-box diagrams, with arbitrary numer- 
ators and integer powers of the propagators were 
developed in [ |7j in the planar case and in [ |^ 
in the non- planar case. In [ the first of the 
two most complicated master integrals involved 
is with all powers of propagators equal to one. 
As a second complicated master integral, the au- 
thors of [ 0] have chosen the diagram with a dot 
on the central line. As was pointed out later [ 
^ , in practical calculations one runs into a linear 
combination of these two master integrals with 
the coefficient 1/e, so that a problem has arisen 
because the calculation of the master integrals in 
one more order in e looked rather nasty. Two 
solutions of this problem immediately appeared. 
In [ |l^ , the authors calculated this very combina- 
tion of the master integrals, while in [|ll[ another 
choice of the master integrals was made: as a 
second complicated master integral, the authors 
have taken the integral which is obtained from the 
first master integral by inserting a specific numer- 
ator. This was a more successful choice because, 
according to the calculational experience, no neg- 
ative powers of e occur as coefficients at these two 
master integrals. 

These analytical algorithms were successfully 
applied to the evaluation of two-loop virtual cor- 
rections to various scattering processes [ |l^ in 
the zero-mass approximation. 

In the case, where one of the external legs is 
on-shell, pi ^ 0, pf = 0, i = 2,3,4, the planar 
double box and one of two possible non-planar 
double-box diagrams with all powers of propa- 
gators equal to one were analytically calculated 
in [ |l^, as functions of the Mandelstam vari- 
ables s and t and the non-zero external momen- 
tum squared pf. Explicit results were expressed 
through (generalized) poly logarithms, up to the 
fourth order, dependent on rational combinations 
of , s and t, and a one- and (in the non-planar 
case) two-dimensional integrals with simple inte- 



grands. To do this, the method based on MB 
integrals mentioned above was applied. These 
and other master planar and non-planar double 
boxes with one leg off-shell were evaluated in [ 
Q with the help of the method of differential 
equations [ The corresponding results are 

expressed through so-called two-dimensional har- 
monic polylogarithms which generalize harmonic 
polylogarithms [ ^ . 

A reduction procedure that provides the pos- 
sibility to express any given Feynman integral to 
the master integrals was also developed in [ Q . It 
is based on the Laporta's observation that, when 
increasing the total dimension of the denominator 
and numerator in Feynman integrals associated 
with the given graph, the total number of IBP and 
Lorentz-invariance equations grows faster than 
the number of independent Feynman integrals. 
These techniques were successfully applied [ |l^ 
to the Feynman integrals with one leg off-shell 
contributing to the process e^e~ — > 3jets. 

For another three-scale calculational problem, 
where all four legs are on-shell and there is a 
non-zero internal mass, a first analytical result 
was obtained in [ |8) for the scalar master double 
box. These and other future similar results will 
be used for calculations connected with Bhabha 



scattering (see also [ 19 for some steps in this 
direction). 

It is believed that sooner or later we shall 
achieve the limit in the process of analytical eval- 
uation of Feynman integrals so that we shall be 
forced to proceed only numerically. (See, e.g., 
[ where this point of view has been empha- 
sized.) However the dramatic recent progress in 
the field of analytical evaluation of Feynman in- 
tegrals shows that we have not yet exhausted our 
abilities in this direction. Indeed, several power- 
ful methods were developed last years. To calcu- 
late the master integrals one can apply the tech- 
nique of MB integration and the method of differ- 
ential equations mentioned above. To construct 
appropriate recursive algorithms one can use re- 
cently developed methods based on shifting di- 
mension [ ^ and differential equations [ Q as 
well as a method based on non-recursive solutions 
of recurrence relations [ |22| . 

One can also hope that new analytical re- 
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suits can be obtained for many other classes of 
Feynman integrals depending on two and three 
scales. In particular, the analytical evaluation 
of any two-loop two-scale Feynman integral with 
two, three and four legs looks quite possible. In 
fact, when going to a higher level of calculational 
complexity, one increases the number of loops, 
legs and independent variables. As the expe- 
rience with the double boxes with one leg off- 
shell has shown, the crucial point is to introduce 
an appropriate class of functions. In this exam- 
ple, these are two-dimensional harmonic polylog- 
arithms which turn out to be adequate functions 
to express result for diagrams of the given fam- 
ily. Presumably, when turning to a situation with 
one more kinematical invariant, e.g. when one 
more leg is off-shell, a natural procedure will be 
to introduce three-dimensional harmonic polylog- 
arithms, etc. 

However, when the number of loops is increased 
one does not need to introduce new functions, so 
that this transition means a 'pure' calculational 
complication. 

2. Triple boxes 

Let us now turn our attention to three-loop on- 
shell massless four-point diagrams and consider 
dimensionally regularized massless on-shell pla- 
nar triple box diagram shown in Fig. 1 . The gen- 
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Figure 1. Planar triple box diagram. 



eral planar triple box Feynman integral without 
numerator takes the form 



T(ai,...,aio;s,i;e) = 
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where s = (pi +^2)^ and t — (^2+^3)^ are Man- 
delstam variables, and fc, I and r are loop mo- 
menta. Usual prescriptions fc^ = fc^ + iO, s ~ 
s -\- lO, etc. are implied. 

By a straightforward generalization of two-loop 
manipulations, one arrives [|2^ at a sevenfold MB 
representation of (|^). In the case of the master 
triple box, = 1, i = 1, 2, . . . , 10, we have 

n(0) 



T^'^^(s,i;e) =r(l,...,l;s,t;e) 
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A"' r(i + w)r(-w) 

s) r(l-2e + w-Z4) 
r(-e + Z2)r(-e + Z3)r(l + w - Zi)V{^z^) 

r(l + Z2 + Z4)r(l+Z3 + Z4) 

^, r(l + e + Z4)r(z2 + Z4)r(z3 + Zi)T{-z^) 

r(i-z5)r(i-z6) 

T{w + Z2 + Z3 + Z4- Z'j)T{-Z2 - zz- Z4) 

T{1 - 2e + Z5 + ze + zy) 
xr(l + e + u> — Z4 — Z5 — — Z7)r(l + Z7) 
xr(— 1 — e — Z5 — Z7)r(— 1 — e — zg — Z7) 
xr(l + Z5 + zg + Z7)r(— e - w — Z2 + Z5 + Z7) 
xr(2 + e + Z5 + Z6 + Z7) 

xr(-e - w - zs + ze + zj) . (2) 

As a first step of three-loop calculations, let us 
consider the evaluation [ p3| of the leading power 
asymptotic behaviour of (j^ in the Regge limit 
t/s — > This calculation demonstrates that a 
three-loop BFKL [ Q analysis is possible. 

One can use the strategy of expansion by re- 
gions [ |2^, |6[ |2^ which shows that in the leading 
power only (Ic-lc-lc) and (2c-2c-2c) regions con- 
tribute, with the leading power 1/t. (See and 
Chapter 8 of [ ^ for definitions of these contribu- 
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tions.) The leading power (2c-2c-2c) contribution 
for the master planar triple box takes the form 



d'^fcd'^/d^r 
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(2pifc + s){2pil + s)(r - l)^P{k - /)2 
1 



(2pir + s) (r + p2 + p)2 



(3) 



where p is such that = t, 2pip = 0, 2p2P = 
—t. The leading power (Ic-lc-lc) contribution 
is obtained due to the symmetry {1^3, 4 ^ 
6, 8 ^ 10}. 

On the other hand, one can organize the cal- 
culational procedure in such a way that the cal- 
culation of the Regge asymptotics is a part of 
the calculation of the unexpanded triple box. 
An analysis of the integrand shows that the key 
gamma functions that are responsible for the 
leading Regge behaviour are r(— e + 2:2,3) and 
r(— 1 — e — Z6,5 — Z7). The standard procedure 
of shifting contours and taking residues can be 
applied. It results again in a sum of MB inte- 
grals where a Laurent expansion of the integrand 
in e is possible. The final result for the Regge 
asymptotics of the planar triple box is [ E3| : 
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where L = Ins/t, 7e is the Euler constant and 
^(z) is the Riemann zeta function. 

For the analytic evaluation of (||), without 
expansion, a similar procedure can be applied. 
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where, up to 1/e^ terms, we have 
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are the generalized polylogarithms [ ^ . 

The coefficient functions , i = 2, . . . , 6, have 
been confirmed [29 by a numerical check with 



the help of numerical integration in the space of 
alpha parameters [ 0. (This algorithm is based 
on a procedure of resolution of singularities in the 
alpha representation that was used in early pa- 
pers on renormalization theory (see, e.g., [ |3^ ) 
and in proofs of similar results for asymptotic ex- 
pansions of Feynman integrals in limits of mo- 
menta and masses (see, e.g., [ 33 and Appendix B 

of 

The present status of the calculation is as fol- 
lows. After integrating over some of the vari- 
ables Zi, ci and cq are expressed through two- 
dimensional MB integrals. In some of them, 
the last integration over a z-variable can be per- 
formed and resulting integrals over w can be eval- 
uated by closing the integration contour to the 
right and taking residues at w = 0, 1, 2, . . .. The 
corresponding results are expressed in terms of 
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harmonic polylogarithms [ |3^ with parameters 
and 1. In the rest of the contributions to Ci and 
Co, subintegrations cannot be done due to Barnes 
lemmas. Hopefully, such two-dimensional MB in- 
tegrals can be analytically evaluated by means of 
'experimental mathematics' similar to used in [ 
p^ , with results presumably expressed in terms 
of harmonic polylogarithms. 

The procedure described above can be applied, 
in a similar way, to the calculation of any massless 
planar on-shell triple box. 
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